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Abstract 

The lowest 37000 eigenvalues of the area operator in loop quantum 
gravity is calculated and studied numerically. We obtain an asymp- 
totical formula for the eigenvalues as a function of their sequential 
number. The multiplicity of the lowest few hundred eigenvalues is 
also determined and the smoothed spectral density is calculated. The 
spectral density is presented for various number of vertices, edges and 
SU(2) representations. A scaling form of spectral density is found, 
being a power law for one vertex, while following an exponential for 
several vertices. The latter case is explained on the basis of the one 
vertex spectral density. 



1 Introduction 

Loop quantum gravity pfl E! amis at quantization of the gravitational field 
without relying on a background metric, while respecting the 3+1 dimen- 
sional nature of spacetime. This approach is based on Ashtekar's new canon- 
ical variables [3] which reduce the constraints to a polynomial form. Note 
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that the configuration space variables are connections rather than compo- 
nents of the metrical tensor. 

Upon quantization, geometrical quantities like area or volume become 
operators acting on the Hilbert space. Area operators are especially interest- 
ing, since their spectral properties rendered possible recently testing the the- 
ory, namely, near black holes correspondence with the Bekenstein-Hawking 
entropyjj] in the quasiclassical limit has been established E]- 

It is a technically important feature of the theory that, by using the so 
called loop representation, the Hilbert space is spanned by the spin network 
states. Such a state is characterized by a graph, whose edges are labelled 
by irreducible SU(2) representations (or, correspondingly, by integer or half- 
integer numbers) and whose vertices lie on a given surface (whose area is to be 
characterized). Vertices are labelled by zero (unit respresentation of SU(2)) if 
gauge invariant states are to be considered (like in this paper). Edges can be 
grouped according to their position with respect to the given surface. Hence, 
an edge can be of type u, d or t, corresponding to up, down and tangential 
orientation. Irreducible representations on edges emanating from the same 
vertex and being of the same type should be multiplied directly and the 
product should be subsequently decomposed to irreducible representations. 
A state is labelled also by one of these irreducible representations (for each 
vertex and each edge type). Thus, to each vertex there belongs a triple 
j u , ja and it- These are not quite independent, as it is required that the 
decomposition of the direct product j u <S> id ® jt should contain the unit 
representation. 

2 Eigenvalues of the area operator 

In the present paper we study the spectral properties of the area operator 
numerically. Actually a lot is known about its eigenvalues and eigenstates 
analytically. [7J Each spin network state is an eigenstate of the area operator. 
The eigenvalues are explicitly expressed as (with a sum over the vertices) 

S = E \ U.A.i,,i + 1) + 2j d>i (j d>i + 1) - j t>i (j t>i + 1) (1) 
1 1=1 

where lp stands for the Planck length and (3 for the Immirzi parameter whose 
value is argued to be jHj . Note that for this latter value the Bekenstein- 
Hawking black hole entropy can be obtained from the theory. 
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Our aim here to study the spectrum numerically, i.e., to make it numeri- 
cally transparent what the formulae yield. As long as just the eigenvalues are 
of interest, it is enough to specify the triple (j u ,i, jd,h jt,i) f° r each vertex, 
the only restriction being that 

\ju,i - jd,i\ < jt,i < ju,i + jd,i (2) 

Hence, j U j and jd,% run over all integer and half-integer number, while 
jt t i runs from \j Uji — to + jd,i in unit steps. In the actual calculation 
a maximal value J for j U)i and jd,% was applied and checked subsequently 
whether the spectrum changed if J was further increased. Moreover, because 
of the symmetric role of j U j and jd,i, it was enough to consider the case when 

Ju,i ^ ]d,i- 

After having the eigenvalues, they have been sorted. The eigenvalues 
versus their sequential number are displayed in Figure If plotting the 




2 

o J 1 1 1 1 1 1 1 

5000 10000 15000 20000 25000 30000 35000 40000 

serial number of eigenvalue 

Figure 1: The lowest 37000 eigenvalues versus their sequential number. 

exponential of the eigenvalues, the curve becomes asymptotically a straight 
line (cf. Fig. |21). Thus, we have asymptotically 
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Figure 2: Fitting the lowest 37000 eigenvalues with a logarithmic function. 



S n = ln(a n + b) . (3) 

with 

a = 0.931888 ±0.0002612 (4) 

b = -2660.66 ± 5.976 (5) 

Certainly, for small sequential numbers n the asymptotical expression 
becomes meaningless. 

The spectrum is discrete, but the mean nearest level spacing decreases for 
increasing eigenvalues. It is also clear, that the spectrum is not bounded from 
above. The largest nearest level spacing occurs between zero (the smallest 
eigenvalue) and the next one, which is 
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3 Density of states of the area operator 



In order to get the density of states, one also needs to know the multiplicities 
of the eigenvalues. This can be calculated for a given number of vertices, 
edges and maximal value J. If the eigenvalue Sj has the multiplicity gj, then 
the density of states may be expressed formally as 



j 

One may smoothen this density by calculating its convolution with some 
function h to get 



The multiplicities were calculated as follows. In the one vertex case the 
irreducible representations on the edges of the same type must be directly 
multiplied and subsequently decomposed to irreducible representations again, 
to get the numbers j u , jd an d it- To get the multiplicity associated with these 
numbers (say, j u ) , one can perform the standard procedure of adding angular 
momenta. This is the following. Suppose that the irreducible representations 
on the type u edges are labelled by ui, u 2 , ... un- I case of one edge, one 
has a multiplicity of 1u\ + 1. For two edges one has to decompose the direct 
product U\ ® u 2 into irreducible representations. The resulting j ranges from 
I ui — u 2 1 to «i + u 2 . The corresponding multiplicity (since a given j appears 
only once) is again 2j + l. For more than two edges one can multiply directly 
the further irreducible representations (one by one) with the decomposition 
obtained in the previous step. Now, however, at a given level one may get 
the same irreducible representation j in several way. Therefore, if we are to 
decompose the product jk-i <8> u k , and a resulting irreducible representation 
is j k (certainly |j fe _i - u k \ < j k < j k -i + u k ), then the multiplicity of j k has 
to be increased by 



P (S) = Y,9AS-Sj) 



(6) 




(7) 




(8) 
, the 
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Note that the multiplicity accociated with the type t edges is simply 2j t + 1, 
because j t has to arise in the decomposition of the product j u ® ja- 

In case of several vertices one proceeds as follows. Multiplicities corre- 
sponding to different vertices must be multiplied together, while the eigenval- 
ues (i.e., yj2j Uti (j Uti + 1) + 2j d>i (j dti + 1) - j t ,i(jt,i + 1)) add U P- If the same 
eigenvalue arises from different graphs, the corresponding multiplicities add 
up. As a consequence, if the density of states is p 1 (S) for one vertex, it will 
be 

J dS 1 p 1 (S - S 1 )p 1 (S 1 ) (10) 

for two vertices, and 

J dSt J dS 2 ..J dS v _ lP 1 (S-S 1 -S 2 -..-S v _ 1 )p 1 {S 1 )p 1 (S 2 )..p 1 {S V - 1 ){ll) 
for V vertices. 

Numerically we obtained that the one vertex density of states p 1 {S) fol- 
lows asymptotically a power law (cf. Fig. |3J), namely, 

p'iS) oc S 6 - 2 , (12) 

roughly independently of J and TV (maximal number of edges). 

For several vertices we obtained that the density of states is asymptoti- 
cally an exponential (cf. Fig. EJ), 

p v {S) ocexp(3.72S) . (13) 

Eq. (J13j) looks a bit contradictory in view of Eq. (jllj) . since a convolution of 
power laws yields power law again. It turns out, however, that one has to 
take into account the discreteness of the spectrum, too. We checked this in 
terms of a simple model of the one vertex density of states. This contained an 
isolated zero eigenvalue, followed by a gap, and then by a continuous power 
law behavior. Taking the convolution of this model density of states, we were 
able to reproduce the approximately exponential behavior of p v {S). 

Note that the proportionality constants did depend on J (maximal index 
of irreducible representations on an edge) and N (maximal number of edges 
emanating from a vertex), but the exponents did not. 
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Figure 3: Density of states for one vertex. 



4 Summary 

We presented a numerical study of the spectrum of the area operator. Espe- 
cially, we calculated the lowest 37000 eigenvalues, and found numerically the 
logarithmic asymtotic formula 0- For the density of states we obtained the 
power law (|T2|) in case of one vertex, and the exponential law (|T3*|) in case 
of several vertices. We explained the connection between these two cases. It 
is remarkable that the spectral properties are rather independent of J and 
N (even of V), exhibiting a kind of scaling behavior. Another interesting 
observation is that in the general (i.e. V > 1) case the density of states 
is exponential, making possible a formulation of thermodynamics, but lead- 
ing to an upper bound for the possible temperatures. Let us note that no 
constraints has been used throughout. 
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Figure 4: Density of states for several vertices. 
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